Simple pendulum dynamics: revisiting the Fourier-based 

approach to the solution 

Riccardo Borghj*] 

Dipartimento di Elettronica Applicata, Universita degli Studi "Roma tre" 
Via della vasca navale 84, 1-00144 Rome, Italy 

Abstract 

The Fourier-based analysis customarily employed to analyze the dynamics of a simple pendulum 
is here revisited to propose an elementary iterative scheme aimed at generating a sequence of 
analytical approximants of the exact law of motion. Each approximant is expressed by a Fourier 
sum whose coefficients are given by suitable linear combinations of Bessel functions, which are 
expected to be more accessible, especially at an undergraduate level, with respect to Jacobian 
elliptic functions. The first three approximants are explicitely obtained and compared with the 
exact solution for typical initial angular positions of the pendulum. In particular, it is shown that, 
at the lowest approximation level, the law of motion of the pendulum turns out to be adequately 
described, up to oscillation amplitudes of 7r/2, by a sinusoidal temporal behaviour with a frequency 
proportional to the square root of the so-called "besinc" function, well known in physical optics. 
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I. INTRODUCTION 



The simple pendulum is the first example, met at an undergraduate level, of a nonlinear 
oscillator whose temporal law of motion does not have a solution expressible in terms of 
elementary functions^, being it customarily given via the use of the so-called Jacobian elliptic 
functions^, not yet available in the "box of tools" 3 of first-year physics students. As a 
consequence, elementary analyses of the pendulum dynamics are generally limited to "small" 
oscillations, for which the differential equation can be linearized to predict an approximately 
harmonic law of motion along the pendulum circular trajectory. Such isochronous solution 
can be intrepreted, through the joint use of a perturbative approach originally proposed by 
Lord Rayleigh^~- and of a Fourier series expansion, as the lowest-order term of a sequence 
of analytical approximants of the exact solution that can be generated iteratively. This has 
been done by Fulchner and Davis 7 , which provided polynomial approximants, with respect 
to the oscillation amplitude, of the Fourier coefficients of the first three elements of the 
sequence. The results obtained in- 7 - have been used in some of the several papers published 
for nearly thirty years to give new and original insights, both theoretical and experimental, 
about the pendulum dynamics^ - — . 

The aim of the present paper is to provide further insights to the Fourier-based treatment 
of the pendulum dynamics. Following the Rayleigh's path, we also propose an iterative 
approach to solve the nonlinear equation in terms of truncated Fourier series arranged in a 
sequence of steps of growing complexity. However, differently from 7 -, the Fourier expanding 
coefficients are not expressed via polynomial approximants but rather obtained, through 
suitable orthogonal expansions, in terms of suitable linear combinations of Bessel functions^. 
A quantitative comparison with the exact solution is then carried out for the first three 
iterations and shows that, at the lowest approximation level the temporal law of motion 
of the pendulum turns out to be adequately described for oscillation amplitudes up to n/2 
(with accuracies better than 0.1 % as far as the pendulum period is concerned) by the so- 
called "besinc" function, well known in optics where it describes mathematically the light 
distribution produced via diffraction by circular holes^ 3 -. 
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II. THEORETICAL ANALYSIS 



A. The Fourier analysis 

We denote by i? = $(t) the angular position of a simple pendulum (with respect to the 
vertical direction) as a function of the time t. Newton's law gives for the well known 
nonlinear differential equation 

i} + col sin ^ = °> (!) 

with uq denoting the natural frequency of the pendulum oscillations and the dots derivation 
with respect to t. We consider symmetric (with respect to the vertical direction $ = 0) 
oscillations having amplitude $o, so that Eq. p]) must be solved together with the initial 
conditions 

tf(0) = #„, 

(2) 

tf(0) = 0. 

Harmonic solutions of this equation are achieved in the limit of small amplitudes (i.e., for 
|$| <C 1) which, via the approximation sin?? ~ transforms Eq. ([T]) into 

$ + u 2 ti = 0, (3) 

whose integral, satisfying initial conditions (J2J), turns out to be 

#(t) = $ cosu t , (4) 

representing isochronous oscillations with period T = 2n/uQ. For amplitude values be- 
yond the above approximation, the pendulum oscillations are no longer isochronous and the 
corresponding period, say T, display the following (exact) dependance on $cr: 



where K(-) denotes the complete elliptic integral of the first kind 2 . Moreover, the exact law 
of motion •& = $(£) is given by 



$0 

$(t) = 2 arcsin 1. sin — sn 
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where the symbol sn(-|-) denotes the Jacobian elliptic function^ (Ch. 22) and u = 2%/T. 
Equation (jSJ) is periodic with respect to t with period T, so that it can be written as a 



Fourier series. Due to the parity of the function sin^ in the r.h.s. of Eq. (JTJ, only odd 
multiples of the fundamental frequency u will appear in the Fourier expansion of d which, 
on taking Eq. (J2J) into account, reads 

oo 
k=0 

where the expanding coefficients c& (k = 1, 3, 5, ...) are functions of the sole initial amplitude 
oscillation $o- Differently from the exact expression given in Eq. ((6]), the Fourier repre- 
sentation in Eq. ([7]) presents a more transparent mathematical structure from which the 
role played by the various expanding coefficients in contributing to give the temporal law of 
motion its temporal "shape" can be inferred. In particular, for "small" initial positions all 
coefficients but the first will turn out to be negligible, being at the same time c\ ~ i9o- All 
Fourier coefficients could be, in principle, obtained starting from the analytical form of 
the solution given in Eq. ([6]) through the following expression: 

c 2fc+ i = | [ cos [{2k + l)ut] dt = 

1 Jo 

(8) 

0o ' 
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(k = 0, 1,2, . . .) which, although cannot be given a closed form, are easily attainable via 
standard numerical packages. In the present paper all numerical calculations have been per- 
formed via the commercial software Mathematica 8.0. The problem of giving analytical 
estimates of the above expanding coefficients has been tackled by Fulchner and Davis 7 , which 
provided Taylor-like polynomial approximants, with respect to i?o> of c\, C3, and C5. 

B. The lowest-order solution 

In the present paper we show further analytical approximants of the coefficients (jHJ) 
that can be obtained by solving Eq. (JTJ) at a different levels of approximations, each of them 
obtainable from the previous one in an iterative manner. To this end we start on substituting 
from Eq. (JTJ) into Eq. (JTJ) so that 

u 2 (ci cos ut + 9c 3 cos3u;t + . . .) = 

(9) 

= Uq sin(ci cosut + C3 cos3ojt + . . .) , 
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and require the expanding coefficients {c2k+i} and the frequency u to be found simply 
on imposing that the periodic functions in both sides of Eq. fl9]) have the same Fourier 
representation, together with the supplementary constraint 



Cl + c 3 + ... = #„, ( 10 ) 

which directly arises from the initial conditions in Eq. (jSJ). We are going to solve Eq. (Q by 
successive approximation, starting from the lowest-order solution which corresponds to the 
choice Ci ~ $0 and c 3 = c 5 = . . . = 0, so that we have 



$0 cos cut ~ Wq sin(-# coswt) , (11) 



where the function in the r.h.s. can be expanded as a Fourier series by using formula (10.12.3) 
of Ref.- to obtain 



sin($ coswi) = 2 (-l) fc J 2k +i(^o) cos [(2k + l)ut] = 

(12) 

= 2Ji(i? Q ) cosut — 2 J 3 ($o) cos3a;t + . . . , 

with J n (-) denoting the nth-order Bessel function of the first kind^. In particular, to find the 
lowest-order estimate of the frequency u it is sufficient to keep in the Fourier expansion (fl2|) 
only the first term which, once substituted into Eq. (fTTT) . gives at once 

uj 2 2Ji(#„) 



> , (13) 

Alternatively we can multiply both sides of Eq. (|lip by cos ut, integrate with respect to the 
dimensionless variable £ = out within the interval [0, 2ir], and recall the integrals 

«2tt 

COS 2 £ d£ = 7T , 







2tt 

cos£ sin(^ cos£)d£ = 27rJ (i? ), 



Before proceeding it worth noting that the function in the r.h.s. of Eq. ( !T3|) is exactly the 
so-called "besinc" function, defined as 

besinc(a) = , (14) 

a 
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FIG. 1: Behaviour, as a function of $o, of the Fourier expanding coefficient c\ numerically evaluated 
from Eq. ([8]) (open circles), together with its lowest-order linear approximation (solid line). For 
completeness the polynomial approximants of c\ provided in 7 is also shown (dotted curve) 

well known in diffraction theory, where describes mathematically the characteristic figure 
produced by light via diffraction through circular holes, the so-called "Airy disk"—. Figured] 
shows the comparison between the exact values (open circles) of the expanding coefficient ci, 
obtained by numerically evaluating the r.h.s of Eq. (|E|) for k = 0, and the lowest-order linear 
approximation (solid line), for ?9 e [0, n]. For completeness the polynomial approximants of 
Ci provided in 7 - is also shown (dotted curve). To give an idea to the reader about the limits 
of the lowest-order approximation of $(t), i.e., 

~ A) coswt, (15) 

with u being given by Eq. (11 3p . Fig. [2] shows, for different values of i?o> a visual comparison 
between the temporal behaviours of the exact (open circles) law of motion in Eq. (|6]), the 
lowest-order estimate in Eq. ( IT5l) (solid curve), and the harmonic solution (dotted curve) 
given in Eq. @. Figure E(a), corresponding to i?o — ?r/3, displays an optimum agreement 
with the exact solution, as well as Fig. Efb) (??o = 77 '/2), for which the agreement, at least 
at a visual level, could be considered acceptable. Clearly the same does not happen for 
figures (c) and (d), corresponding to $0 = 37r/4 and 97r/10, respectively, for which a deeper 
investigation is required. 
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FIG. 2: Visual comparison between the temporal behaviours of the exact (open circles) law of mo- 
tion, the lowest-order estimate in Eq. f)15[) (solid curve), and the harmonic solution in Eq. @ (dot- 
ted curve), for 0„ = vr/3 (a), vr/2 (b), 3vr/4 (c), and 9vr/10 (d). 



C. Higher-order solutions 

The derivation of the lowest-order solution can be further refined in the following way: 
the first-order estimate of i?(t) will obtained by truncating the Fourier series in Eq. ([7]) up to 
k = 1, i.e., by including now the first two harmonics. To find the corresponding coefficients, 
say {c^c'a}, and also the new estimate of the frequency we can substitute the lowest-order 
solution given in Eq. (fT5l) in the r.h.s. of Eq. (Q and to expand it, thanks to the Fourier 
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series expansion in Eq. ([12]) . up to the term containing cos3o;t to obtain 



to 



— (c[ COS tot + 9 C3 cos 3ut) 



to.. 



(16) 



~ 2 Ji(-^o) costot — 2 J 3 (i? ) cos3wt, 
where also the frequency to has to be re-calculated by using again Eq. (TTDT) written as 



ci + 4 ~ i? 



• 



(17) 



so that, after some algebra, we obtain 



L l — 



c 3 — 



to 
tol 



where the correction term rj is given by 



1 + 77 



?9o 

1 + i 
rj 



2 MA 



(18) 



(1 + V) 



1 ^3(^0) 



^"9» (19) 
Similarly as done in Fig. [H Fig. [3] shows the behaviour of the expanding coefficients 
given in Eq. (fT8|) as functions of the initial amplitude #o- Moreover, Fig. [3] shows the same 
comparison reported in Fig. [2j which now displays an agreement with the exact solution 
that is really good for the cases depicted in figures (a), (b), and (c), but it is still far from 
being satisfactory for the case of figure (d). The reason is due to the fact that for the latter 
case the pendulum starts from an initial position close to the vertical, which leads to a law 
of motion with a "flat" shape in the space (t, $) (this is more evident for values $0 which 
are closer to it) and, accordingly, it would require a greater number of harmonics to achieve 
a better representation. We prefer not to detail here the whole development of the second- 
order solution, which is confined in the El for interested readers, but to limit ourselves to 
note that it will be obtained by keeping in the Fourier series (I7j) the first three terms whose 
expanding coefficient estimates, say {c'(,c^,c^} are shown, as functions of $0, in Fig. [5]^a), 
(b), and (c), respectively. Figure |6] shows the results obtained by using the second-order 
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FIG. 3: The same as in Fig. [1] but for the first-order solution, (a): Behaviour of the coefficient c^. 
(b): Behaviour of the coefficient c' s . 

solution, which includes the first three harmonics in the Fourier expansion of Eq. ([7]). In 
particular figure (a) corresponds to the choice $0 — 97r/10, in order for a visual comparison 
with Fig. 3(d) to be allowed. It is seen that the second-order solution displays an excellent 
agreement with the exact solution considerably better than that corresponding to the first- 
order one. However, for initial positions closer to the vertical also the second-order solution 
reveals to be inadequate, as shown in figures (b) and (c), so that one should pass to the 
third-order solution, and so on. 

To conclude the present section we wish to compare the estimates of the pendulum 
period, obtained at the above described three subsequent approximation levels, with the 
exact expression. This is shown in Fig. [7(a), where the exact behavior of T/Tq as a function 
of 6q is shown (open circles) together with the estimate provided by the lowest-order (dashed 
curve), the first-order (dotted curve), and the second-order (solid curve) approximants of 
the law of motion. To give a more quantitative information, in Fig. UJb) the corresponding 
values of the relative error (in percentage) are also plotted. 

III. CONCLUSIONS 



Nonlinear physics is a fascinating world where it is not always possible to establish general 
rules similarly to what happens for most linear problems, that can often be solved through 
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FIG. 4: The same as in Fig. [2] but with the first-order solution in place of the lowest-order one. 
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FIG. 5: The same as in Fig. [IJ but for the second-order solution, (a): Behaviour of the coefficient 
c'{. (b): Behaviour of the coefficient C3. (c): Behaviour of the coefficient C5. 



"closed- form approaches." The study of the simple pendulum dynamics is not an exception 
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FIG. 6: Comparison between the temporal behaviours of the exact (open circles) law of motion 
and the second-order estimate (solid curve) for $o = 97r/10 (a), 957r/100 (b), and 997r/100 (c). 
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FIG. 7: (a): exact behavior of T /Tq as a function of 9q (open circles), together with that corre- 
sponding to the lowest-order (dashed curve), the first-order (dotted curve), and the second-order 
(solid curve) solutions, (b): corresponding values of the relative error (in percentage). 

and it is really difficult to find elementary ways to approach it and to offer its solution to 
undergraduates. An approach to the solution of the simple pendulum problem which should 
appear to be "natural" is to exploit the periodicity of the temporal law of motion by repre- 
senting it as a Fourier series, with unkown expanding coefficients that should be found via 
the constraints imposed by the nonlinear differential equation and the appropriate initial 
conditions. In the present paper such an approach has been used to build up analytical 
approximants of the pendulum solution, as well as of the period, in terms of suitable super- 
positions of Bessel functions of the first kind. In particular, at the lowest approximation level 
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the pendulum period turns out to be well described (with accuracies better than 0.2 %for 
oscillation amplitudes up to 7r/2) by the so-called "besinc" function, well known in physical 
optics. We also believe that, from a mere didactical point of view, such a connection could 
be exploited to convince even an audience of skeptic students about the fact that different 
phenomena in physics are often, depending on the approximation level, quantitatively de- 
scribed by equations having similar mathematical aspects 34 . Accordingly, methods aimed 
at finding their solutions must necessarily share a common signature. This, in particular, 
is true as far as the key role played by some classes of so-called "special functions" (and 
the Bessel family is one of the most important) is concerned, a fact that makes them of 
fundamental importance in applications to several different branches of physics 35 . 
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Appendix A: Derivation of the second-order solution 

Again we have to start from Eq. (Q written in the following form: 

u 2 

— (c'[ cos cut + 9c 3 ' cos3cut + 25 c'l cos5cut)) ~ 
cu 

(Al) 

~ sin(c' 1 cos cut + c 3 cos3cut) , 
where {c'^Cg} denotes the set of coefficients found at the first level and given in Eqs. f fT8|) 
and (fT9|) . while the triplet {c'{, c 3 ', c'l} denotes the set of coefficients that we have to include 
in the representation of $(t) at the second level of expansion. To find them consider the 
Fourier series expansion of the function sin i9(t) at the first level, namely 

sm(c[ cos cut + c' 2 cos3cut) = 

(A2) 

= C\ cos cos cut + C 3 cos cos 3cut + C 5 cos cos 5cut + . . . , 
which will be truncated up to the first three terms before inserting into Eq. flAlj) with the 
expanding coefficients {Ci,C 3 ,C 5 } being given by 

1 f 2n 

C2/C+1 = — cos [(2k + l)x] sin(c' 1 cosx + c' 3 cos3x). (A3) 
Wo 
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Although the above integrals cannot be expressed exactly in closed form, analytical estimates 
of them can be obtained by approximating the integrand as follows: 

sin(c^ cos a; + c' 3 cos3x) ~ sinfc^ cosx) + c' 3 cos 3a; cos(c' 1 cosx) , (A4) 

which, once substituted into Eq. (IA3j) . after some algebra gives 

1 f 2lT 

C 2 k+i — — cos [(2/c + l)x] sin(c 1 cosx) 
Wo 

(A5) 

c' f 2lT 

+ — cos [(2k + l)x] cos 3x cos(c[ cosx) , 

7T Jo 

and, on taking formulas (10.12.2) and (10.12.3) of^ into account, 

(-l) h C 2k+1 ^ 2J 2k+1 (c[) + c' 3 [J 2k+4 (c[) - J 2 k-2{c' 1 )\. (A6) 



Finally, on substituting from Eq. (1A2|) into Eq. ( 1A1[) and on comparing both sides of the 
resulting equation term by term, the following system is obtained: 

i^-WfW' k = 0,1,2, (A7) 
which, together with the boundary condition 

4 + 4 + 4 = $o , (A8) 

and on taking Eq. (IA6[) into account, allows both the frequency uj and the triplet {c", c' 3 ', C5} 
to be straightforwardly found. 

Of course the above described procedure can be further iterated to retrieve, in principle, 
trigonometric approximants of the pendulum law of motion up to arbitrarily high accuracies, 
that is left to the interested reader. 



13 



* Electronic address: borghi@uniroma3.it 

1 A. Sommerfeld, Lectures on Theoretical Physics. I. Mechanics (Academic Press, 1970) 

2 NIST Digital Library of Mathematical Functions. |http://dlmf.nist.gov/ , Release 1.0.5 of 2012- 
10-01. 

3 R. P. Feynamn and R. Leighton, Surely you're joking, Mr. Feynman! (Norton & Co., NY, 
1985). 

4 Lord Rayleigh, The Theory of Sound (Cambridge University Press, Cambridge, 1899). 

5 N. Minorsky, Introduction to Nonlinear Mechanics (Edwards, Ann Arbor, MI, 1947). 

6 H. Lamb, Dynamics (Cambridge University Press, Cambridge, 1961). 

7 L. P. Fulchner abd B. F. Davis, "Theoretical and experimental study of the motion of the simple 
pendulum," Am. J. Phys. 44, 51 - 55 (1975). 

8 R. Simon and R. P. Riesz, "Large amplitude pendulum: a Fourier analysis," Am. J. Phys. 47, 
898 - 899 (1979). 

9 S. C. Zilio, "Measurement and analysis of large-angle pendulum motion," Am. J. Phys. 50, 450 

- 452 (1981). 

10 W. P. Ganley, "Simple pendulum approximation," Am. J. Phys. 53, 73 - 76 (1985). 

11 L. H. Cadwell and E. R. Boyco, "Linearization of the simple pendulum," Am. J. Phys. 59, 979 

- 81 (1991). 

12 M. I. Molina, "Simple linearizations of the simple pendulum for any amplitude," Phys. Teach. 
35, 489 - 90 (1997). 

13 R. B. Kidd and S. L. Fogg, "A simple formula for the large-angle pendulum period," Phys. 
Teach. 40, 81 - 83 (2002). 

14 L. E. Millet, "The large-angle pendulum period," Phys. Teach. 41, 162 - 163 (2003). 

15 R. R. Parwani, "An approximate expression for the large angle period of a simple pendulum," 
Eur. J. Phys. 25, 37 - 39 (2004). 

16 G. E. Hite, "Approximations for the period of a simple pendulum," Phys. Teach. 43, 290 - 292 
(2005). 

17 A. Belendez, A. Hernandez, T. Belendez, A. Marquez, and C. Neipp, "Analytical approximations 
for the period of a nonlinear pendulum," Eur. J. Phys. 27, 539- 551 (2006). 



14 



18 F. M. S. Lima and P. Arun, "An accurate formula for the period of a simple pendulum oscillating 
beyond the small angle regime," Am. J. Phys. 74, 892 (2006). 

19 A. Belendez, A. Hernandez, T. Belendez, C. Neipp, and A. Marquez, "Application of the ho- 
motopy perturbation method to the nonlinear pendulum," Eur. J. Phys. 28, 93 - 104 (2007). 

20 P. Amore, M. C. Valdovinos, G. Orneles, and S. Z. Barajas, "The nonlinear pendulum: formulae 
for the large amplitude period," Rev. Mex. Fis. E 53, 106 - 111 (2007). 

21 S. Gil, A. E. Legarreta, D. E. Di Gregorio, "Measuring anharmonicity in a large amplitude 
pendulum," Am. J. Phys. 76, 843 - 847 (2008). 

22 C. G. Carvalhaes and P. Suppes, "Approximations for the period of the simple pendulum based 
on the arithmetic-geometric mean," Am. J. Phys. 76, 1150 - 1154 (2008). 

23 F. M. S. Lima, "Simple log formulae for pendulum motion valid for any amplitude," Eur. J. 
Phys. 29, 1091 - 1098 (2008). 

24 A. Belendez, J. J. Rodes, T. Belendez, and A. Hernandez, "Approximation for a large-angle 
simple pendulum period," Eur. J. Phys. 30, L25 - L28 (2009). 

25 Yuan Qing-Xin and Ding Pei, "Comment on 'Approximation for a large-angle simple pendulum 
period'," Eur. J. Phys. 30, L79 - L82 (2009). 

26 A. Belendez, J. Frances, M. Ortuho, S. Gallego, and Jose G. Bernabeu, "Higher accurate 
approximate solutions for the simple pendulum in terms of elementary functions," Eur. J. Phys. 
31, L65 - L70 (2010). 

27 M. Turkyilmazoglu, "Improvements in the approximate formulae for the period of the simple 
pendulum," Eur. J. Phys. 31, 1007 (2010). 

28 M. I. Qureshi, M. Rafat, and S. Ismail Azad, "The exact equation of motion of a simple 
pendulum of arbitrary amplitude: a hyper geometric approach," Eur. J. Phys. 31, 1485 (2010). 

29 K. Ochs "A comprehensive analytical solution of the nonlinear pendulum," Eur. J. Phys. 32, 
479 (2011). 

30 Augusto Belendez, Enrique Arribas, Andres Marquez, Manuel Ortuho, and Sergi Gallego, "Ap- 
proximate expressions for the period of a simple pendulum using a Taylor series expansion," 
Eur. J. Phys. 32, 1303 (2011). 

31 T. G. Douvropoulos, "Simple analytic formula for the period of the nonlinear pendulum via the 
Struve function: connection to acoustical impedance matching," Eur. J. Phys. 33, 207 (2012). 

32 Xue De-Sheng, Zhou Zhao, and Gao Mei-Zhen, "Accurate Period Approximation for Any Simple 

15 



Pendulum Amplitude," Chinese Phys. Lett. 29, 044601 (2012). 

M. Born and E. Wolf, Principles of Optics (Cambridge University Press, Cambridge, Mas- 
sachusetts, 1999). 

See for instance the beautiful discussion at Chapter 12 of the second volume of the Feynamn 
Lectures on Physics about the analogies in physics. 

M. V. Berry, "Why are special function special?," Phys. Today 4, 11-12 (2001). 



16 



